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Global existence and asymptotic behavior of solntions to a nonlocal 

Fisher-KPP type problem 


Shell Bian* Li Chen'^ Evangelos A. Latos^ 


Abstract 

In this work, we consider a nonlocal Fisher-KPP reaction-diffusion problem with Neumann 
boundary condition and nonnegative initial data in a bounded domain in IR”(n > 1), with 
reaction term m“(1 — m{t)), where m{t) is the total mass at time t. When a > 1 and the initial 
mass is greater than or equal to one, the problem has a unique nonnegative classical solution. 
While if the initial mass is less than one, then the problem admits a unique global solution for 
n = 1, 2 with any l<a<2orn>3 with any 1 < a < 1 -|- 2/n. Moreover, the asymptotic 
convergence to the solution of the heat equation is proved. Finally, some numerical simulations 
in dimensions n = 1, 2 are exhibited. Especially, for a > 2 and the initial mass is less than one, 
our numerical results show that the solution exists globally in time and the mass tends to one 
as time goes to infinity. 


1 Introduction 


In this work we consider the following nonlocal initial boundary value problem. 


ut — Alt = it“ 

_G 

1 

1—1 

X G Q,t > 0, 

(la) 

Vit • = 0, 


X G dQ, 

(lb) 

u{x, 0) = uo{x 

) >0, 

X G 17, 

(Ic) 

where u is the density, 17 is a 

smooth bounded domain in W^, 

n > 1, a > 1 and 

is the outer 


unit normal vector on 50. Without loss of generality, throughout this paper we assume |0| = 1 
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(otherwise, rescale the problem by |n|), let m{t) = f^u(x,t)dx and mo = m(0). A damping term 
with (T > 0 can also be included to get ut — Au + au = rt“(l — m{t)). In this case, similar results 
to this paper can also be obtained. For simplicity, we assume that a = 0. 

In the 1930s, Fisher m and Kolmogorov, Petrovskii, Piskunov |23j in population dynamics and 
Zeldovich, Frank-Kamenetskii m in combustion theory started to study problems with this kind 
of reaction terms. Actually, they introduced the scalar reaction-diffusion equation ^ -\-F{u), 

and studied the existence, stability and speed of propagation. In the theory of population dynamics, 
the function F is considered as the rate of the reproduction of the population. It is usually of the 
form 

F{u) = /3n“(l — u) — ju. 

From the above model two cases emerge depending on the values of a. 

In the case of a = 1, the reproduction rate is proportional to the density u of the population 
and to available resources {1 — u). The last term, —'yu, describes the mortality of the population. 

The case a = 2, which is the motivation for our work, considers the addition of sexual repro¬ 
duction to the model with the reproduction rate proportional to the square of the density, see [l2] . 
For more information on reaction-diffusion waves in biology, we refer to the review paper of Volpert 
and Petrovskii m- 

Next we pass to the relation between the local and the non-local consumption of resources. In 
the local reaction-diffusion problem 

du d'^u ^ . 

where u is the population density, describes the random displacement of the individuals of 
this population and the reaction term represents their reproduction and mortality. Moreover, the 
reaction term consists of the reproduction term which is represented by the population density to 
a power, ti", multiplied with the term (1 — u) which stands for the local consumption of available 
resources. 

The nonlocal version of the above problem is 

du d'^u / \ 

^+ “ y -y)'^{v^t)dyj -1'^^ (3) 

where /S ,7 > 0 and = 1. It can be seen as the case where the individual, located 

at a certain point, can consume resources in some area around that point. (j){x — y) represents 
the probability density function that describes the distribution of individuals around their average 
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positions and it depends on the distance from the average point x to the actual point y. One can 
easily verify that if (/> is a Dirac 5-function, then the nonlocal problem reduces to Q. 

In the current paper we will study problems with reaction terms similar to the above nonlocal 
reaction terms. There are some already known results on the reaction-diffusion equation with a 
nonlocal term, 

ut = Au + F{t,u,I{u)), I{u) = / u{y,t)dy, 

Jo, 

in a bounded domain D. However, compared to the local version, the results for the nonlocal 
reaction terms of Fisher-KPP type are relatively limited. Here we list some of the known recent 
results. 

Anguiano, Kloeden and Lorenz considered F = f{u)I{u‘^){l — I{u^)) and proved the existence 
of a global attractor [T]. Wang and Wo |44) proved the convergence to a stationary solution with 
F = — I{vF), m > 1. For F = -|- 6/(e'>'“), Pao [33] studied the existence or nonexistence 

of stationary solutions. Liu, Chen and Lu |28j proved also the blow-up of solutions for a similar 
equation, see also [Hj . Rouchon obtained global estimates of solutions in m- For more information 
on nonlocal KPP-Fisher type problems, we refer to a recent book by Volpert [40] . 

Nonlocal Fisher-KPP type reaction terms can describe also Darwinian evolution of a structured 
population density or the behavior of cancer cells with therapy as well as polychemotherapy and 
chemotherapy, we refer the interested reader to the models found in |29l ISO] El] . 

Bebernes and Bressan [6] (see also Bebernes 0, Pao|33]) considered the equation with reaction 
term 

F{t, u,I{u)) = f{t, u{t, x)) -I- [ g{t, u{t, y))dy, t > 0, x E D. 

Jn 

They considered the case when f{t,u) = e“, g{t,u) = ke^ {k > 0), for which the above problem 
represents an ignition model for a compressible reactive gas, and proved that solutions blow-up. 
Later, Wang and Wang |45] considered a power-like nonlinearity, i.e. 

F{t,u,I{u))= / vF{t,y)dy — ku^{t,x), t > 0, x E D, 

Jn 

with p,q > 1, and proved the blow-up of the solutions. 

Budd, Dold and Stuart Hu and Yin [22] considered a similar to the above problem in the 
case p = 2 and general p respectively, 

F{t,u,I{u)) = uP - j^uP{t,y)dy, t > 0, x E D. 

With this typical structure, the energy of the solutions is conserved (under Neumann boundary 
conditions). For this kind of nonlocal problems it is known [l5] that there is no comparison principle 
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and they are the closest models to the ones we are considering in this work. For a general study on 
nonlocal problems, we refer to Quittner’s and Souplet’s book m as well as the paper by Souplet 

m- 

1.1 Preliminary discussion 

In this article, we will focus on 0 which has a reaction term of the type 

F{t,u, I{u)) = ^ — J u{x,t)dx 

for a > 1. One additional fact that makes this problem more difficult to handle is the lack of 
comparison principle as one can see for example from |45j . 

For nonnegative u, formally by integrating Q over 11, we get, 

m{t) = (1 — m{t)) / u°‘dx, 

Jn 

where m{t) = u{x,t)dx is the total mass at time t. 

If we start at time to such that 1 — m{to) < 0, which means that m(to) > 1, we can see that 
m'{t) is negative and therefore m{t) decreases in time. In this case it is natural to expect the global 
existence of solutions. 

On the other hand, if we start at time to such that 1 — m{to) > 0, we can see that m{t) increases 
in time. However if (1 — m{t)) remains positive, the equation has a similar structure to the heat 
equation with a power-like reaction term for which we know that the problem might have no global 
solution for super-critical exponent a < 1 + 2/n (see for example [U [5l |T8l [IHl |20l [23] ). In this paper 
we give a negative answer to this observation. Our main results are the following two theorems: 

Theorem 1. Let n > 1, a > 1 and f^uo(x)dx = mo > 0. Assume uo is nonnegative and 
uo G L^{Q) for any 1 <k < oo. Then for mo < 1 with a satisfying 

1 <a <l + 2/n, n>3, (4) 

l<a<2, n = l,2, (5) 

or mo > 1 with arbitrary a > 1, problem & has a unique nonnegative classical solution. Moreover, 
the following a priori estimates hold true. That’s for mo < I, 

\W{-,t)\\Lk(^^) <C+ C t~^ foranyt>0. (6) 
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For mo > 1 


n > 3, 
n = 1, 2. 


(7) 

( 8 ) 


Here C denote different constants depending on mo,k,a, but not depending on ||wo||Lfc(n)- 

Theorem 2. Letu{x,t) be the unique nonnegative classical solution obtained from Theorem^ v be 
the solution to the heat equation with Neumann boundary condition and initial data f^VQ{x)dx = 
mo, then as t ^ oo, 

||?x(-,t) - v{-,t) - {mo - 1)||l2(j^) < Cie~'^^\ (9) 

where Ci,C 2 are constants depending on the initial mass mo and ||ifo||L2“(n)• 

This paper is organized as follows. In Section 2, we firstly present the dynamics of the mass. 
The global existence of the solutions and thus the proof the Theorem are shown in Section 3. 
Section 4 is devoted to the proof of Theorem Section 5 shows the numerical simulations of the 
problem which give the motivation for our future study for the case a > 1 + 2/n. Finally, Section 
6 concludes the main work of our paper, some open questions for problem Q are also addressed. 

2 Dynamics of the total mass 

The evolution of the total mass plays the key role in our proof for the global existence of the classical 
solution. Therefore, we firstly give the evolution of mass u{t)dx in time. 

Lemma 3. For mo > 0, the mass f^u(t)dx = m{t) satisfies 

min{l,mo} < m{t) < max{l,mo}. (10) 

Furthermore, we have the following decay estimates 

|1 - rn{t)\ < |1 - (11) 

Proof. We return to the original problem 0 and integrate it over 0 to get: 

m'ft) = {1 — m{t)) I u'^dx. (12) 

Jn 

There are two possibilities depending on the initial mass. 
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• If we start at time to where 1 — m{to) > 0, we can see that m![t) is positive and therefore m 
increases in time. Moreover, with the use of Jensen’s inequality we get by using m{t) > mo, 

fn'{t) > (1 — m{t))m°‘(t) > (1 — m{t))mQ. 

By solving this inequality we get a lower bound on the speed with which m{t) increases to 1 

i.e. 

m{t) > 1 — — mo). 


• If 1 — m(to) < 0, then m decreases in time. By monotonicity, we get m{t) < mo and again by 
Jensen’s inequality, 

m {t) < (1 — m(t))m“(t) < 1 — m{t), 

then we get that 

m{t) < 1 + (mo — l)e“*. 


By putting together the above two cases we have the expected results. □ 


3 Global Existence 


This section mainly focuses on the global existence of the classical solution to 0. We will use the 
following ODE inequality from [9] , which was also used in [10] . 

Lemma 4. Assume y{t) >t) is a function for t > 0 satisfying 

y'{t) <a- j3y{tY 


for a > 1, a > 0, ,0 > 0, then y{t) has the following hyper-contractive property 

1 


y{t)<{a/l5f/- + 

Furthermore, ify{0) is bounded, then 


fi{a — l)t 


for any t > 0. 


y{t) < max (y{0), . 


(13) 


(14) 


The proof of global existence heavily depends on the following two a priori estimates, Proposition 
and Proposition and then we will use the compactness arguments to close the proof. Due to 
the preliminary discussion, we will divide the a priori estimates into two cases mo < 1 and mg > 1. 
Firstly, we focus on mo < 1. 
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Proposition 5. Let n > 1 and mo <1. If a satisfies 


1 < a < 1 + 2/n, n > 3, 

1 < a < 2, n = 1, 2, 

then for any I < k < oo, the nonnegative solution of 0 satisfies 

k-l 

\W{'fi)\\L'‘{n)^C{mo,k,a) + C{mo,k,a)t for any t > 0. (15) 

Moreover, if uo{x) G Lf{Q), then 

ll«(•>^)llLfe(^) < max|||'Uo(a:)||Lfe(Q),C(mo, A:,a)| , (16) 

and for any 0 < T < oo 


Vnt G (0,T;L^(fl)) . 

Proof. Since mo < 1, by lemma one has mo < m{t) < 1. Using ku^~^ as a test function for 
equation Q and integrating it by parts 


— / u'^dx = - 

dtJn 

dtJn k 


4(A; - 1) 


k 


f \'Vu2\'^dx + k f ^dx(l— f udxj, 

Jo. Ju \ ) 


f \Vu^\‘^dx + km{t) f ^dx = k 

Jo, Jo, 


u 


k-\-a—l 


dx. 


(17) 


Choosing l<A:'<fc + a — 1, combining Holder’s inequality and the Sobolev embedding theorem 
one has 


in 


/-Ini/ / A fc 2(/e+Q: —1) /, y\k‘2ik+oc — l) 


(l__^) 2(fc+a-l) 




2(k + a-l) 


< 

ti2 

k 

U2 



Lp{n) 



2k' 

L^{n) 


<C{k) 

<C{k) 


Vut 




A 

k 

1-A 

k 

A 

k 

L2(Q) 

M2 

2k' + 

Lnr(n) 

M2 

L2(Q) 

M2 


2A(fc+t:i —1) 
k 

L2(n) 


k 

rt2 


2(1-A)(fc + Q;-1) 
k 

2k' 

L^{n) 


+ 


k 

U2 


1-A 

2k' 

Ln^{n) 

2\(k + a-l) 
k 

L‘^{n) 


2{k + a-l) 
k 


k 

U2 


2(l-A)(fc+c«-l) ' 
k 

2k' 

Lur(n) 


(18) 


where A is the exponent from Holder’s inequality, i.e. 

k k 


A= £(0,1), 

2k' p 


(19) 
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and p satisfies 


P=^2^ ^>3, 

2 (fc+a-l) ^ n = 2, 


p = oo, 


n = 1. 


Now we will divide the analysis into three cases n > 3, n = 2 and n = 1. 


For n > 3, p = and then 


n—2 


kn kn 


\ _ 2(/c+a 1) ^ /'n 1 ^ 

- hr, , ^ I'J’ 


kn _|_ 1 n 

2k' ^ 2 


with k > max I ^ Taking k' > simple computations arrive at 


n\n \ kn I (a—l)n 

l\\k + 0 — 1) pp - p- n 


k 


kn _|_ 1 _ n 

2k' ^ 2 


< 2 . 


To sum up, for k' > max | 2 '^^ ■> ’ thanks to the Young’s inequality, from (18) one has 

k-l 


/n k 


Vwt 


+ Cik) 


tt2 


2 

L2(n) 

. 2A(fc+i: 

+ C{k) 

r-l) 

k 

L2(n) 


k 

tt2 


k 

U2 


(l_^) 2(fc+a-l) 1 


2k' 

L^{n) 

2(l-\)(k + a-l) 
k 

2k' 

L^(0) 


A(fc + Q! —1) 
^ k 


Letting 


r = (1 - A) 


2{k + a-l) 


k i_ 

recalling m{t) > mo, together 0 with ( |22[ ) we arrive at 

d 

In 


1 

A(/c+a—1) ’ 


^ / vl^dx + kniQ I —^||Vrta ||^2(n) 

dt Jfi Jq k 

On the other hand, using Holder’s inequality with 1 < k' < k + a — 1 we have 


l-V 


where 


{k + a — l){k' — 1) 


A:'(A: + a — 2) 


. , (A: + a - 1)(A: - 1) 

kikJ-2) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 

(26) 


(27) 










































Hence 


kr 

I 2 




kr 


krO 




( 28 ) 


Taking (24)-(28) into account we obtain that 


d 


ii„.i|fe+o-l , Iiv7„,# I|2 


— I^u dx + kmo\\u\\J+c-i(Q) + 


(O) 


k 






(29) 


Here 


h = (i-Ap+«-i)—4^. 

k 

k ot — 1 k “ 1 “ (y. — 2 

Recalling the definition of 6, r], A, direct computations show that Xrj + (1 — \)9 < 1. For 


1 ^ Oi ^ 1 T —, 
n 


(30) 


one can derive that 


krO 

< k -\- a — 1. 


(31) 


Next for n = 2, < p < oo, by proceeding the similar arguments to the case n > 3 from 


(21) to (29), we obtain that for 


1 < a < 2 - 

P 


(32) 


(31) holds true. When n = 1, p = oo, then for 1 < a < 2, (31) also holds true. 

Therefore, combining the three cases n > 3, n = 2 and n = 1, using Young’s inequality we 


obtain from (29) that 


j^u’^dx + kmo\\u\\lt+Ji(Q) + 


< . +C{mo,k). 


4 ' 4 

In addition. Holder’s inequality yields that 


' k 
iLfc+“-i(n) 


L2{n) 


(33) 


u 


LH^)J 


^ ||^,||fc+o-l 


a — 1 
I k-1 




(34) 
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Then using lemma we solve the ODE inequality 


dt Jn 


a-l 

k-1 


9IU; II ^ 


1+1 


u^dx I < C{mo, k) 


(35) 


to obtain that for any 1 < k < oo, 

I k 


\u\ 


Lk{p,) — C{mQ, k, a) + 


C{mo, k, a) 


k-1 
ct — l 


for any t > 0. 


(36) 


Furthermore, if uo{x) G for any 1 < A; < oo, then taking y{t) = u^dx in lemmaj^one has 

||ii||Lfe(Q) < max|||uo(ai)||Lfc(f7),C'("io, A:,a)| . (37) 


Now we integrate (33) from 0 to T in time, then we can obtain that for any T > 0 

rT r r 


/ u^{T)dx+ / / 

/o Jo Jn 


Vn2 


dxdt + / / ^dxdt < / UqJx + C{mo,k)T, 

Jo Jn Jn 


from which we derive 


Vut G (0, T; L^(D)) for any T > 0. 


This completes the proof. □ 

For mo > 1, owing to lemma we know m{t) > 1 for any t > 0, thus we have the following 
result 

Proposition 6. Let n > 1 and a > 1. Assume uq G L)^(D) and f^uo(x)dx = mo > 1, Then for 
any 1 < k < oo, the nonnegative solution of satisfies that for any t > 0 


lik^n) — k) + 

lifc(Q) < C{mo, k) + 


C{mo,k) 


n(fc —1) 


, n > 3, 


C{mo,k) 


1 k-1 


, n = l,2. 


Moreover, if uq G L^{Ll), then 


roo p 

7o Jn 


Vu2 


dxdt < 


f u^dx. 

Jn 


(38) 

(39) 

(40) 

k 


Proof. Recalling lemma we know that if mo > 1, then m{t) > 1 for any t > 0. Hence the L 
estimates 0 can be reduced to 


dt/n k 


/ |Vu2 \'^dx < 0. 

Jn 


(41) 
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For n > 1, using Holder’s inequality one has that for any 1 < k < oo, the following estimate holds 

\k ^ \\„.\\ke 


\U 


= llip(n)ll'“llLi(Q)^ 


(42) 


where 6 = i- and 
fe—- 

p 


V = — 
■e n-2’ 


n >3, 


2 < p < oo, n = 2, 
p = oo, n = 1. 


(43) 


Thanks to the Sobolev embedding theorem and Young’s inequality, from (42) one has 
\k 


U\ 


L'=(n) 


< C(n) (^||Vu 2 + ||'U2 |Il2(o) 

2 


\u\ 


k{i-e) 

I ^ 
iLi(n) 


< C(mo,n)||Vu2||^2(j^) + C{mo,n)\\u\\’lk^^^ 

< C{mo,n)\\Vu^\\l 2 ^^) + ^ (hllii(n))" +C{mo,n,k). (44) 


Plugging the above estimates into (41) yields that 


^ J u^dx + C{mo,n,k) u’^dx ^ + C{mo,n,k)\\Vu^\\j^ 2 (^^') < C{mo,k,n), (45) 


solving the ODE inequality we have that for any t > 0 


u 


\u 


lifc(n) < C{mo,k,n) + 
< C{mo,k,n) + 


C{mQ, k, n) 
t 

C(mo, k, n) 
t 


n(fc —1) 


1 k-l 


n> 3, 
n = 1,2. 


Moreover, if uq G T*(H), then (41) directly yields 


l^llL'=(r2) < ll'*^o||L'=(r2)- 


Next integrating (41) from 0 to oo in time we obtain that 

roo r 

Jo Jn 


Vu2 


dxdt < / u^dx. 


(46) 

(47) 

(48) 

(49) 


This closes the proof. □ 


Remark 7. In fact, (38) and (39) also hold true for heat equation, and the uniform boundedness 
in time of the norm depends only on the initial mass, not depends on the initial norm. 
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Now we have obtained the necessary a priori estimates to complete the proof of Theorem It 
can be proved by standard methods and for the convenience of the reader we mention the key steps 
in the following. First of all, from the above estimates, we can take k = 2 and fc = 2a in Proposition 
[^andj^to get the estimates for ||Vm||j;^2(^2(o,t)) \\ut\\L^(^H-'^{o,T}) for T > Q. By Aubin-Lions 

lemma nulla] , we have the strong compactness of u in so that the nonlinear terms can be 
handled. Therefore, the global existence of weak solutions (in the sense of distributions) can be 
obtained by standard compactness argument. Secondly, from the estimates of the weak solution in 
Propositionand Proposition]^ the nonlinear term u“(l — m{t)) G L^([0,T] x 12), V/c > 1 for any 
r > 0. The solution is a strong solution from classical parabolic theory, [261127j . By Sobolev 
embedding, we can bootstrap it to get that classical solution. In the end, the uniqueness can be 
obtained directly from comparison principle 123123 since ^(1 ~ rn{t)) is bounded from below. 
□ 


4 The long time behavior of solutions 

As we can see from the above arguments, equation Q has a unique classical solution. In this 
section, we will detect the long time behavior of the global solution. 

Theorem 8. Assume uq G for any 1 < k < oo. Let u he the classical solution to problem 

0 and V be the solution to the heat equation with Neumann boundary condition and initial data vq 
such that f^vo(x)dx = mo, Then as t ^ oo 

l|w(-T) - - (mo - 1)||l2(j^) < Cie~^^\ 

where the constants Ci, C 2 depend on mo and ||uo||L 2 a(o). 

Proof. The difference between the two equations is 


(n — v)t + A{u — v) = u“(l — m{t)). 


Let u{t) = f^u{x,t)dx and v{t) = f^v(x,t)dx. By (12), we have ut = m'{t) = (1 — m{t)) f^u°‘dx. 
vt(t) = 0 because of v{t) = vo- Therefore, 

{u — v)t — {u — v)t + A(u — v) = n“(l — m(t)) — (1 — m(t)) / n“. 

Jn 

The standard estimate shows that. 


blT f \{u—v)—{u—v)\‘^dx+f \V{u—v)fdx = {l—m{t)) 
2 dt Jq Jq 


u — u 

Jn 


u 


— v) — {u — vf) 


dx. 
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By taking k = 2a m Proposition and Proposition we get 

f \{u — v) — {u— v)\‘^dx + [ \V{u — v)fdx 
2 dt Jq Jq 

< ^|1 —m(t)| / \{u — v) — (u— v)\'^dx + C\1 — m{t)\, 

2 Jn 

where C depend on mg, a and ||ifo||L2“(n)- Applying Poincare inequality and lemmawe get 

hX [ — v) — {u— v)\'^dx + C{^}) f \{u — v) — {u— v)\‘^dx 

2 dt Jq Jq 


-ct 


< Ce / \{u — v) — {u — v)\^dx + Ce 
Jn 

From the above ODE we get the following estimate, 

f \{u — v) — {u— v)\‘^dx < 

Jn 

where Ci,C 2 are constants depending on mo and ||uo||l 2 c«(q). Thus completes Theorem□ 

5 Numerical Results 

For 1 < a < 2 (n = 1, 2), if mo < 1, it has been shown in the previous sections that the solution 
will exist globally without any restriction on the initial data. An interesting question is whether 
the solution also exists globally for a > 2. In the following, we will give a complete numerical study 
for n = 1, 2 with any a > 1. 

5.1 Numerical scheme 

For the numerical simulation, we consider the 2-dimensional equation in D = [0, b]'^ with any 6 > 0 

ut = Au -j- udxdy) , (x, y) G [0, b] x [0, 6], t> 0, 

u{t = 0) = uo{x, y) > 0, (x, y) G [0, b] x [0, b], 

= ^ib,y,t) = 0, 

. = ^{x,b,t) = 0. 

Denote 


(50) 


f{u) = tt" ( 1 — I udxdy 

In 


(51) 
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Here an alternating direction implicit difference scheme is applied to construct the numerical com¬ 
putations. 

Let h be the space step and r be the time step, N = h/h + 1 \s the number of the discrete points, 
T is the final time and K = T/t + 1. Denote 


x{i) = {i — 1) * h, i = 1,2,3, ■ ■ ■ , N — 1, N, 
y{j) = (j - 1) * j = 1, 2, 3, • • • , - 1, N, 

t{k) = {k — 1) * T, k = 1,2,-■ ■ ,T/t + 1, 

= [{x{i),y{j))\l < i,j < 

The discrete solution at each time is presented as a matrix ■ G , where 

Ui,j = u{x{i),y{j),^, l<i,j<N, 

= u(^-,-,t{k)^, k = l,2,---,K. 


Next we introduce some notations: 


x2„, _ 2r(ij + Ui-ij 

O^Uij — ^2 : 


We construct the ADI scheme with operator splitting method and using Taylor expansion lin¬ 
earizes f{u) to approximate this semi-linear equation as follows 


- u^. 


= K + s: 


% ^ % 


{« +1) Si ( 


- Uij 


^^+1 _ yk 

+ f{0+ '' o ,2<i,j<N-l, 


ujj = Uo (x{i),y{j)'^, l<i,j <N. 


(52) can be rearranged into 


(52) 




1 a2 
1 2^y 


= h{u. 


u / ’ 


h () = 




rT-2 _2 

' -2 . ' 




(53) 
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u(x.y.t) when t=2s 



(a) u{x, y) at t=0.5s 


(b) u{x, y) at t=2s 




(c) u{x, y) at t=7s 


(d) max of u{x, y) with time evolution 


Figure 1: u{x,y) with time evolution, initial mass mo < 1 
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(a) mass of u{x, y) with time evolution 
Figure 2: Mass u with time evolution, initial mass mo < 1 


(53) is equivalent to the following ADI scheme 


= 2<i<N-l, j = 2,3,--- ,1V-1, A: = l,2,-.. ,K. (54) 


4u2j — U3J _ 4tt7V-l,j — UN-2,j 
Ul,j = - 5 -, UN,j = -^- 


, j = 2,3,--- ,iV-l, 


1- = 2<j<N-l, i = 2,3,--- ,N-1, k = l,2,--- ,K-1, 


4y,^+l _ 7/^+^ 
..k+l _ ^^*.2 ^i,3 

^i,i - 3 


fc+1 _ ^i,N-2 


, z = 2,3,--- ,iV-l. 


5.2 Numerical examples 

Let D = [0,1]^. In the numerical study, we will consider three cases: 

1. n = 2,1 < a < 2, we choose mo < 1 and uq G L^{Vt,). 

2. n = 2,a > 2, we choose mo < 1 and uq G L"(D), J^UqcIx is large enough. 

3. n = l,a = 2, we choose mo < 1 and uq G L“(D), f^u^dx is large enough. 

Case 1: n = 2, we choose a = 3/2 and the initial data 

uo{x, y) = (—2x^ + 3x^ + 0.5) (—2y^ + 3?/^) , / / uo{x, y)dxdy = 0.5 < 1. 


(55) 

(56) 

(57) 


(58) 

Fig. [T] shows the evolution of the solutions with time. We notice that the solution converges to 
1/|D|, where |D| is the area of the 2-n domain. Fig. [^shows its corresponding total mass with time 
evolution. Eventually the mass converges to 1. 
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(a) u with time evolution 



(b) max of u with time evolution 

Figure 3: u{x, y) with time evolution, initial mass mo > 1 
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(a) u{x,y)dxdy with time evolution 

Figure 4: Mass of u{x, y) with time evolution, initial mass tuq > 1 


For the initial mass is greater than 1, we choose 


uo{x, y) = {-2x^ + 3x^ + l) (-2/ + + l) , 


// uo{x,y)dxdy 


2.25 > 1. 


(59) 


j Jn 

the results are shown in Fig. where the solution converges to l/|f2|. The evolution of mass is 
shown in Fig. we observe that it finally converges to 1. 

Case 2: n = 2, a = 3 and mo < 1, the initial data is chosen to be a characteristic function 


u = 


0.3 < x < 0.3 + 5h, 0.3 < y < 0.3 + 5/i, 
0, other, 


where h = 0.01. Simple computations deduce u^dx = mg = 0.625 < 1 and Ugdx = 3.9 x 10^. 
Fig. [5] shows the evolution of the solutions with time. We can notice that the mass tends to one 
very quickly and u{x,y) goes to l/|n| with time evolution. 

Case 3: n = 1, a = 2 and mo < 1, the initial data is chosen to be 


10, 0.3 < X < 0.35, 
0, other. 


(60) 


the initial mass mg = 0.5. It can be observed in Fig. [^that u{x) converges to the constant 1/|12|, 
Fig. [7] shows the mass tends to one finally. 
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(a) initial u{x, y) 


(b) u{x,y) at t — 10 '‘s 




(c) mass of u{x, y) with time evolution (d) max of u{x, y) with time evolution 

Figure 5: n = 2, u{x, y) with time evolution, niQ < 1 
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(a) u{x) with time evolution 


(b) max of u{x) with time evolution 


Figure 6: n = 1, u{x) with time evolution, mo < 1 



(a) u(x)dx with time evolution 

Figure 7: n = 1, mass of u(x) with time evolution, initial mass mo < 1 
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6 Conclusions 


This paper concerns the nonlocal Fisher KPP problem Q with the reaction term’s power a. For 
a > 1, the global existence of a classical solution to 0 is analyzed. When the initial mass is less 
than one, 1 < a < 2 for n = l,2orl<a<l + 2/n for n > 3, there exists a global unique 
nonnegative classical solution. When the initial mass is greater than or equal to one, the Fisher 
KPP problem admits a unique classical solution for any a > 1. Our numerical simulations show 
that when the initial mass is less than one and a > 2 for n = 1,2, the unique nonnegative classical 
solution will exist globally. Therefore, for n = 1,2 with a > 2 or n > 3 with a > 2/n, our 
conjecture is that the problem Q also admits a global unique nonnegative classical solution. This 
is a challenging problem which we will study in the future. 
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